finish
/clear
/filn, E331 
/title, Buckling analysis of short cylinder
! ********* 趨*********
D=1000                         !Ͳо
R=D/2                          !Ͳа뾶
L=5000                         !Ͳ㳤
t=10                           !Ͳ
kmodel=1                       !ģͿأ1ȫģͣ21/2ģ
nhoop=20                       !Բܷ߶εʷ
*if,kmodel,eq,1,then           !ģ趨Ͳʷ
    naxial=40
*else
    naxial=20
*endif
p=1                            !λѹ
EXX1=2e5                       !ϵĵģ
NUXY1=0.3                      !ϵĲɱ

!****************ǰ***************************
/PREP7  
ET,1,SHELL63                   !趨Ԫ
R,1,t,t,t,t, , ,               !ǵԪʵ
MP,EX,1,EXX1                   !ϵĵģ
MP,PRXY,1,NUXY1                !ϵĲɱ

!****************ģ***************************
k,1,0,0,0                      !Ͳ߹ؼ
*if,kmodel,eq,1,then           !ģ趨ģͲĳ
    L=L	
*else
    L=L/2                      !1/2ģͲ峤ȡ
*endif
k,2,0,0,L                      !Ͳһؼ
l,1,2                          !
LESIZE,all, , ,naxial,, , , ,1 !趨Ͳʷ
circle,1,r                     !ͲԲ
lsel,s,loc,z,0                 !ѡԲ߸߶
LESIZE,all, , ,nhoop,, , , ,1  !趨Բʷ
ADRAG,all, , , , , ,1          !Բ
AMESH,all                      !Բʷ

!********************趨Լ*************************
*if,kmodel,eq,1,then           !ģ趨Լ
    nsel,s,loc,z,0             !ȫģͣԼĳһλ
    nsel,r,loc,x,R	    
    D,all,uz,                  !ȫģͣԼĳһλ
    nsel,s,loc,z,0             !ѡһڵ
    nsel,a,loc,z,L             !ѡһڵ
*else
    nsel,s,loc,z,0             !ڰģͣڶԳʩӶԳԼ
    DL,all,,symm
    nsel,s,loc,z,L             !ѡһڵ
*endif
csys,1                         !
NROTAT,all                     !תѡڵĽڵϵ
D,all,UY                       !ʩӻλԼԱ֤ý汣ΪԲ
fini
!**************** ȡ***************************
/SOLU
ANTYPE,STATIC                  !̬
PSTRES,ON                      !ԤӦӰ
SFA,all,2,PRES,p               !Բǵʩѹ
Alls                           !ȫѡ
SOLVE                          !
FINI                           !˳
!**************** ȡ***************************
/SOLU                          !½
ANTYPE,BUCKLE                  !ָΪ
BUCOPT,LANB,1                  !Block Lanczosн㣬ȡһģ̬
MXPAND,1                       !趨ģ̬չΪ1
SOLVE                          !
FINISH
!*******************************************
/POST1
SET,FIRST
PLDISP,1
